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METHOD AND APPARATUS FOR REDUCING COUPLING BETWEEN 

SIGNALS 

Background of the Invention 

Field of the Invention 

The present invention relates to the field of signal processing, and, more particularly, 
relates to multi-channel demodulators for demodulating mixed signals, such as, for example, 
signals generated in a pulse oximetry system. 
Description of the Related Art 

In many multi-channel measurement and communication systems, crosstalk between 
channels and corruption of data within the channels are significant problems. Such problems 
can arise from variations in manufacturing tolerances, movement, propagation delays, phase 
shifts, temperature effects, degradation of components due to age or other factors, noise, etc. 

A pulse oximetry system is one example of a system where the above-referenced 
problems are found. In a pulse oximetry system, blood oxygen saturation is determined by 
transmitting pulses of electromagnetic energy through a portion of a subject having blood 
flowing therein (e.g., through a finger, through an ear lobe, or other portion of the body 
where blood flows close to the skin). The pulses of electromagnetic energy comprise 
periodic pulses of red light having wavelengths of approximately 660 nanometers, for 
example, and periodic pulses of infrared light having wavelengths of approximately 905 
nanometers. 

After propagating through the portion of the subject, the red pulses and the infrared 
pulses are detected by a detector which is responsive to light at both wavelengths and which 
generates an electrical signal that has a relationship to the intensity of the electromagnetic 
energy incident on the detector. The detector output is a two-channel signal having a first 
signal component corresponding to the detected red pulses and a second signal component 
corresponding to the detected infrared pulses. 

The two-channel signal is demodulated to recover separate signals corresponding to 
the first signal component and the second signal component. However, prior art 
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demodulators are not sufficiently accurate enough to completely separate the two signal 
components in all cases. Thus, it is not uncommon for the first demodulator output signal 
(corresponding to the first signal component) to contain residual components of the second 
signal and vice versa. This crosstalk between the first and second signal components reduces 
the accuracy of the recovered first and second signals. In multi-channel systems with more 
than two channels, crosstalk can occur between all of the channels, again reducing accuracy. 

Summary of the Invention 
The present invention solves these and other problems by separating a combined 
multi-channel signal into a plurality of output signals in a manner that reduces crosstalk and 
other contamination in the plurality of output signals. In one embodiment, the separator 
includes a multi-channel demodulator that is first configured using nominal values for the 
various components in the signal path. In one embodiment, the multi-channel demodulator 
is further configured using data obtained from calibration measurements. In one 
embodiment, the calibration measurements are made during an initialization period. In one 
embodiment, the calibration measurements are made frequently, continuously, or at selected 
intervals. In one embodiment, calibrations are performed on at least one of, initialization, on 
command, on attachment of a new sensor, continuously, and/or interspersed with 
measurements. In one embodiment of a system for measuring one or more blood 
constituents, the calibration measurements are made when the system detects that a patient 
has been connected to the system. In one embodiment, the multi-channel demodulator is 
further configured at regular intervals by re-running the calibration measurements. In one 
embodiment, the multi-channel demodulator comprises an optimizing demodulator. In one 
embodiment, crosstalk in the multi-channel demodulator is reduced by computing an 
amplitude and/or phase adjustment of one or more demodulation signals that are provided 
respectively to one or more mixers. 

In one embodiment, an apparatus for measuring blood oxygenation in a subject 
includes a first signal source which applies a first input signal during a first time interval. 
A second signal source applies a second input signal during a second time interval. A 
detector detects a first parametric signal responsive to the first input signal passing 
through a portion of the subject having blood therein and detects a second parametric 
signal responsive to the second input signal passing through the portion of the subject. 



The detector generates a detector output signal responsive to the first and second 
parametric signals. A signal processor receives the detector output signal and demodulates 
the detector output signal by applying a first demodulation signal to a signal responsive to 
the detector output signal to generate a first demodulator output signal and applying a 
second demodulation signal to the signal responsive to the detector output signal to generate 
a second demodulator output signal. In one embodiment, the first demodulation signal has at 
least one component comprising a first frequency, a first phase, and a first amplitude; and 
the second demodulation signal has at least one component comprising a second frequency, 
a second phase, and a second amplitude. In one embodiment, the first phase and the second 
phase are chosen to reduce crosstalk from the first parametric signal to the second 
demodulator output signal and to reduce crosstalk from the second parametric signal to the 
first demodulator output signal. In one embodiment, the first amplitude and the second 
amplitude are chosen to reduce crosstalk from the first parametric signal to the second 
demodulator output signal and to reduce crosstalk from the second parametric signal to the 
first demodulator output signal. In one embodiment, at least one of the first amplitude, the 
first phase, the second amplitude, and the second phase are chosen to reduce crosstalk from 
the first parametric signal to the second demodulator output signal and to reduce crosstalk 
from the second parametric signal to the first demodulator output signal. 

In one embodiment, at least one of the first amplitude, the first phase, the second 
amplitude, and the second phase is determined by turning off one of the first and second 
signal sources and measuring the crosstalk between one of the parametric signals and the 
non-corresponding output signal. 

One embodiment includes a method of reducing crosstalk between two signals 
generated by applying a first pulse and a second pulse to measure a parameter. The first 
pulse and the second pulse are applied periodically at a repetition rate defining a period. The 
first pulse is generated during a first interval in each period and the second pulse is 
generated during a second interval in each period. In one embodiment, the second interval is 
spaced apart from the first interval. In one embodiment, the second interval overlaps at least 
a portion of the first interval. The first and second pulses produce first and second 
parametric signals responsive to the parameter. The first and second parametric signals are 
received by a detector which outputs a composite signal responsive to the first and second 



parametric signals. The method includes applying a first demodulation signal to the 
composite signal to generate a first demodulated output signal. The first demodulation signal 
includes at least one component having at least a first amplitude and a first phase. The 
method further includes applying a second demodulation signal to the composite signal to 
generate a second demodulated output signal. The second demodulation signal includes at 
least one component having at least a second amplitude and a second phase. The method 
further includes lowpass filtering the first demodulated output signal to generate a first 
recovered output signal responsive to the first parametric signal, and lowpass filtering the 
second demodulated output signal to generate a second recovered output signal responsive 
to the second parametric signal. The method also includes choosing at least one of the first 
phase, the first amplitude, the second phase, and the second amplitude to reduce crosstalk 
components in the first recovered output signal and the second recovered output signal. In 
one embodiment, the method also includes choosing the first phase and/or the second phase 
to reduce crosstalk components in the first recovered output signal and the second recovered 
output signal. 

In one embodiment, the first phase and the second phase are chosen by applying a 
first light pulse pattern during a first time period and measuring the first recovered output 
during the first time period as a first calibration output, and measuring the second recovered 
output during the first time period as a second calibration output. The method includes 
applying a second light pulse pattern during a second time period and measuring the first 
recovered output during the first time period as a third calibration output and measuring the 
second recovered output during the second time period as a fourth calibration output. The 
method further includes computing the first phase and the second phase from at least the 
first calibration output, the second calibration output, the third calibration output, and the 
fourth calibration output. 

In one embodiment the first phase is computed from a ratio of the first calibration 
output and the second calibration output. 

In one embodiment, the first demodulation signal includes a sum of a first 
demodulation component having a first amplitude and a second demodulation component 
having a second amplitude. The second demodulation component is in quadrature with the 
first demodulation component and the act of choosing the first phase involves choosing the 



first amplitude and the second amplitude. In one embodiment the quadrature components 
are sinusoidal and cosinusoidal. 

In one embodiment, the first demodulation signal includes a sum of a sinusoidal 
component having a first amplitude and a cosinusoidal component having a second 
amplitude. The first amplitude and the second amplitude are chosen by a least squares 
minimization of an error corresponding to the crosstalk. In one embodiment, the error is 
integrated over a time period corresponding to an integer number of cycles of the sinusoidal 
component. 

In one embodiment, a first demodulation signal is applied to a composite signal 
having first and second coefficients to generate a first demodulated signal. The first 
demodulation signal includes a first component having a first amplitude and a second 
component having a second amplitude. The first and second components being in 
quadrature. The second amplitude has a predetermined relationship to the first amplitude. 
The predetermined relationship is selected to cause the first demodulated signal to have 
lower frequency components that include a primary component corresponding primarily to 
the first desired component and a residual component corresponding to the second 
component. The first demodulated signal is lowpass filtered to generate a first output signal. 
At least one of the first amplitude and the second amplitude are adjusted to reduce the 
residual component with respect to the primary component. 

In one embodiment, a pulse oximetry system includes a modulation signal generator. 
The modulation signal generator generates a first modulation signal including a first pulse at 
a repetition frequency having a first duty cycle. The modulation signal generator generates a 
second modulation signal including a second pulse which also repeats at the repetition 
frequency and having a second duty cycle. The second pulse can be non-overlapping with 
respect to the first pulse, or the second pulse can partially or completely overlap the first 
pulse. The first and second pulses include a plurality of components wherein a first 
component has a frequency corresponding to the repetition frequency and a second 
component has a second frequency corresponding to twice the first frequency. A first 
transmitter emits electromagnetic energy at a first wavelength in response to the first pulse. 
A second transmitter emits electromagnetic energy at a second wavelength in response to 
the second pulse. A detector receives electromagnetic energy at the first and second 



wavelengths after passing through a portion of a subject. The detector generates a detector 
output signal responsive to the received electromagnetic energy. The detector output signal 
includes a signal component responsive to attenuation of the electromagnetic energy at the 
first wavelength and a signal component responsive to attenuation of the electromagnetic 
energy at the second wavelength. A first demodulator multiplies the detector signal by a first 
demodulation signal and generates a first demodulated output signal. A second demodulator 
multiplies the detector signal by a second demodulation signal and generates a second 
demodulated output signal. A configuration module configures the first demodulation signal 
and the second demodulation signal to substantially separate the first demodulator output 
and the second demodulator output. 

In one embodiment, the configuration module selects a phase relationship between 
the first demodulation signal and the second demodulation signal. 

In one embodiment, the configuration module configures the first demodulation 
signal and the second demodulation signal using, at least in part, data obtained during a 
calibration period. In one embodiment, the calibration data includes first and second 
calibration data corresponding to the first and second demodulated output signals during a 
first time period, and third and fourth calibration data corresponding to the first and second 
demodulated output signals during a second time period. In one embodiment, the second 
transmitter is turned off during the first time period, and the first transmitter is turned off 
during the second time period. 

In one embodiment, the configuration module configures the first demodulation 
signal and the second demodulation signal by adjusting initial parameters that define the first 
demodulation signal and the second demodulation signal. The configuration module adjusts 
the initial parameters using, at least in part, the calibration data obtained during a calibration 
period. 



Brief Description of the Drawings 
The present invention will be described below in connection with the accompanying 

figures. 

Figure 1 is a block diagram of a multi-channel processing system that uses feedback 
from one or more outputs to configure the operation of a signal separator that separates a 
composite signal into a plurality of output signals. 

Figure 2 is a block diagram of a two-channel signal processing system to determine 
blood oxygen saturation in a subject, wherein illumination is provided by back-to-back 
Light-Emitting Diodes (LEDs). 

Figure 3 is a block diagram of a multi-channel signal processing system to determine 
blood constituents (e.g., oxygen saturation) in a subject, wherein illumination is provided by 
N diodes or illumination sources. 

Figure 4 is a block diagram of a specific embodiment of the multi-channel 
processing system of Figure 1. 

Detailed Description of the Preferred Embodiment 
Figure 1 shows a topology of a multi-channel measurement or communication 
system 100. The system 100 has a signal combiner 103 for combining one or more input 
signals S\...Sn into a composite signal and a signal separator 104 for separating the 
composite signal into one or more output signals S x . . . S M The output signals S, . . . S M can 
include estimates of the input signals S\ . . .£ N . The input signals S\ . . .SW are corrupted by pre- 
combination distortion 101-102 respectively, and, optionally, by combination distortion in 
the signal combiner 103. The combiner 103 combines the N input signals into a composite 
signal (or composite signals). The combiner 401 can combine signals by addition, 
subtraction, multiplication, division, modulation, non-linear processes, linear processes, 
estimation, combinations thereof, etc. The composite signal is provided through a 
communication channel to the separator 104. The composite signal is distorted by 
communication channel distortion 110. The separator 104 separates the composite signal 
into M output signals, where M can be less than N 9 greater than N, or equal to N. In one 
embodiment, the separator also provides one or more additional output signals n Q ...n K 



corresponding to estimates of other signals, such as, for example, noise signals, error 
signals, etc. 

Due to errors in the system 100, the output signals S { . . . S M are typically not exact 

copies of the input signals, but rather are estimates of the input signals. The accuracy of 
these estimates is a measure of system performance. The pre-combination distortion 101- 
102, the combiner distortion, and/or the channel distortion 110 tend to introduce crosstalk 
between the channels and thereby corrupt the output signals. The pre-combination distortion 
101-102, combiner distortion, and the channel distortion 110 can be caused by variations in 
manufacturing tolerances, delay, movement, temperature effects, degradation of components 
due to age or other factors, noise, etc. 

A module 105 is provided to configure the separator 104 to improve the quality of 
the separation function and thereby improve the quality of the output signals. One or more 
of the output signals from the separator are provided to the module 105 to provide feedback 
regarding the quality of the output signals and/or feedback regarding the operation of the 
separator 104. The module 105 uses feedback from one or more of the output signals 
S, . . . S M (and, optionally, the output signals n Q ...n K ) to monitor the quality of the separation 
function and to provide control information to control the operation of the separator. In one 
embodiment, the module 105 is configured by using configuration data obtained from the 
combiner 103. Such configuration data can be obtained by calibration procedures that test 
the operation of the combiner 103 before or during system use. 

In one embodiment, the module 105 configures demodulators in the signal separator 
104 using, at least in part, calibration data obtained during a calibration period. For example, 
in one embodiment involving a two channel system, the calibration data includes first and 
second calibration data corresponding to the first and second output signals during a first 
time period, and third and fourth calibration data corresponding to the first and second 
demodulated output signals during a second time period. In one embodiment, the second 
transmitter is turned off during the first time period, and the first transmitter is turned off 
during the second time period. In one embodiment, the module 105 configures the first 
demodulation signal and the second demodulation signal by adjusting initial parameters that 
define the first demodulation signal and the second demodulation signal. The configuration 
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module adjusts the initial parameters using, at least in part, the calibration data obtained 
during a calibration period. 

Figure 2 is a block diagram of a two-channel signal processing system 200 that fits 
the general topology shown in Figure 1 . The system 200 is configured to determine one or 
more blood constituents in a subject, such as, for example, a human subject. In the example 
presented, the measurements are performed on a portion of the subject, such as a finger 202 
illustrated in Figure 2 A. An LED modulation circuit 204 drives a pair of back-to-back light 
emitting diodes (LEDs) 206, 208 by applying a periodic signal to the two light emitting 
diodes 206, 208. Light from the diodes 206, 208 passes through the finger 202 and is 
detected by a detector 250. An output from the detector 250 is provided to a signal 
processing block 270. A control output from the signal processing block 270 is provided to 

the LED modulation circuit 204. The signal processing block 270 also provides outputs 5, 
and S 2 corresponding to the light detected from the diodes 206, 208, and, optionally, output 
signals n 0 ...fi K corresponding to estimates of noise or other signals. 

In one embodiment, the LED 206 is selected to emit electromagnetic energy in the 
red visible light range, and has a wavelength of, for example, approximately 660 
nanometers. The LED 208 is selected to emit electromagnetic energy in the infrared range, 
and has a wavelength of, for example, approximately 905 nanometers. The LED 
modulation circuit 204 supplies current to activate the LEDs 206 and 208. Each LED is 
activated for a time period x which can be different for the different LEDs. The pulses from 
the LEDs 206 and 208 repeat with a periodicity T. 

Figure 3 is a block diagram of a multi-channel signal processing system 300 that 
also fits the topology shown in Figure 1. Like the system 200, the system 300 is configured 
to determine blood oxygen saturation or other blood constituents in a subject, such as, for 
example, a human subject. In Figure 3, an LED modulation circuit 314 drives N diodes, 
where N is two or greater, thus allowing greater flexibility than the system 200. In Figure 3, 
the diodes 206 and 208 are shown, along with an N'th diode 309, it being understood that 
the diode 309 is omitted if N=2. The LED modulation circuit 314 is configured to allow the 
diodes 206, 208, and 309 to be driven independently such that the diodes can be driven at 
separate times or in overlapping time periods if desired. Light from the diodes 206, 208, 309 



passes through the finger 202 and is detected by the detector 250. The output from the 
detector 250 is provided to a signal processing block 371. A control output from the signal 
processing block 371 is provided to the LED modulation circuit 314. The signal processing 
block 371 also provides outputs Si through Sm, where M is greater than or equal to one, but 
need not be equal to N. 

Figure 4 shows one embodiment of an adjustable multi-channel 
modulator/demodulator system 400. The system 400 is based on the topology of the system 
100 and can be used in a system for measuring blood constituents (e.g., pulse oximetry, 
carboxyhemoglobin, etc.) as shown in Figures 2 and 3. In the system 400, an input S\(t) and 
a modulation input M\(t) are provided to a first modulator 491. A signal input Sn(0 and a 
modulation input M N (/) are provided to an N lh modulator 402. 

The photodetector 250 is modeled as an adder 405. The outputs of the modulators 
491 and 402 are added together in the adder 405, in the presence of noise n(t) to generate a 
composite signal M(t) where: 

S(t) = S, {t)M x (t)+... + S N (t)M N (t) + n(0 (1) 

The S(t) signal output of the adder 405 (i.e., the output of the detector 250) is applied 
to the input of a signal-processing block 410. Within the signal-processing block 410, the 
signal S(t) is passed through an amplifier 497 and through an analog bandpass filter 498. 
The analog bandpass filter 498 provides anti-aliasing and removal of low frequency noise 
and DC. The desired signal components in the signals S,{t) are frequency shifted by the 
operation of the modulation signals M t {t) and are passed by the analog bandpass filter 498. 

The output of the analog bandpass filter 498 is sampled by an analog-to-digital 
converter 499 and converted therein to digital signals and provided to an input of an optional 
decimation block 420. 

The filtered (and, optionally, decimated) signal S(t) is sampled to produce a 
sampled-data signal S(k) that is provided to the first input of the first mixer 424, to the first 
input of the N 'th mixer 412, and to the first input of a noise channel mixer 413. A first 
demodulating signal D\(k) is provided to a second input of a first mixer 424 from a signal 
generator 431 . The N th demodulating signal Dj^k) is provided to an N* 1 mixer 412 from an 
output of a signal generator 432. The noise demodulating signal D 0 (k) is provided to the 
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noise channel mixer 413 from an output of a signal generator 441. A control input to each 
of the signal generators 431, 432, and 441 is provided by the output of the adjuster algorithm 
450. In yet another embodiment, the adjuster algorithm 450 may also be controlled by other 
signal processing elements downstream of the signal processor 400. 

The outputs of the mixers 413, 424, and 412 are provided as respective inputs to 
decimation blocks 440, 430, and 434 respectively. Each of the decimation blocks 440, 430, 
and 434 has a control input provided by the output of the adjuster algorithm block 450. The 
output of the decimation block 440 is an estimate of the signal n{t) and it is provided to an 
input of the adjuster algorithm block 450. In an alternate embodiment, the signal estimates 

S-Xk) are also provided to the adjuster algorithm block 450. 

An output of the decimator 430 is a signal S x (k) , which, as discussed above, is an 
estimate of the signal S\(k) (where S\(k) corresponds to a sampled-data representation of 
S\(t)). Likewise, the output of the decimation block 434 is an estimate of the signal S/^f). 
As shown above, the selection of the demodulating signals D t {t) for i = 0...N in accordance 
with the present invention substantially reduces or eliminates the effects of noise in the 

output signals S ( (k) and n(k), and also substantially reduces or eliminates crosstalk between 
the signals. 

When the system 400 is used in connection with a blood constituent measurement 
system as shown in Figures 2 and 3 the red LED 206 provides a light intensity represented 
as Ird, and the infrared LED 208 provides a light intensity represented as Iir. The effects of 
turning the LEDs 206, 208 on and off on a periodic bases are modeled by the first multiplier 
or modulator 290 which applies a first modulation signal M\(i) to the red light intensity to 
generate a modulated red signal Irdmod(0 and by a second multiplier the modulator 292 
which applies a second modulation signal Mi(i) to the infrared light intensity to generate a 
modulated infrared signal Iirmod(0- The modulated light red signal and the modulated 
infrared signal are applied to the finger 202, or other body portion, as described above. The 
blood in the finger 202 has a volume and scattering components, which vary throughout 
each cardiac cycle. The blood carries oxygen and other materials therein. The oxygen 
content is a function of both the blood volume and the concentration of the oxygen in the 
blood volume. The concentration of the oxygen in the blood volume is generally measured 
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as blood oxygen saturation for reasons which are described in full in U.S. Patent Nos. 
5,482,036 and 5,490,505, both of which are hereby incorporated by reference in their 
entirety. As further described in the two referenced patents, the blood oxygen saturation is 
determined by comparing the relative absorption of the red light and the infrared light in the 
finger 202. The comparison is complicated by the noise caused by movement, ambient 
light, light scattering, and other factors. The signals S\(t) and S 2 (t) represent the effect of the 
time-varying volume and scattering components of the blood in the finger 202 on the red 
light and the infrared light, respectively, passing through the finger 202 from the LEDs 206, 
208 to the detector 250. 

As shown in Figure 4, a set of N+l signals Si[k] i = 1..JV, and n(k) are sampled at a 
desired sample rate. The signals are combined according to the formula: 

S(k) = Mi(k)Si(k) + . . .+M N (k)S N (k) + n(k) (2) 

In one embodiment, each of the decimators 420, 440, 430, and 434 includes a digital 
lowpass filter and a sample rate compressor. In one embodiment, the characteristics of the 
digital lowpass filters (e.g., the number of filter coefficients and values of the filter 
coefficients) and the sample rate compression factor of each decimator are fixed. In one 
embodiment, the characteristics of the digital lowpass filters (e.g., the number of filter 
coefficients or values of the filter coefficients) and the sample rate compression factor of 
each decimator are provided by the adjustment algorithm 450. The signal generators 431, 
432 and 441 generate the demodulation sequences for the demodulators 424, 412, and 413 
respectively. The demodulation sequences produced by the signal generators 431, 432 and 
441 are controlled by the adjuster algorithm 450. 

In one embodiment, the adjuster algorithm 450 adjusts the pre-demodulation 
decimation rate R\ (in the demodulator 420), and the post-demodulation decimation rate R 2 
(in the demodulators 430, 434 and 440) according to the noise in the noise estimate n(k) and 

(optionally) according to the signals S i (k) . The product R\R2 is the total decimation rate 

from the signal S(k) at the output of the A/D converter 499 to the signals S t (k) at the output 

of the signal processing block 400. The adjuster algorithm may adjust R\ and R2 such that 
the product R\R 2 varies, or the adjuster algorithm may adjust R\ and R 2 such that the product 
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R\R 2 is substantially constant. Typically, the adjuster algorithm will keep the R\R 2 product 
constant so that the signal processing blocks downstream of the signal processor 400 will 
operate at a substantially constant sample rate. 

In one embodiment, the adjuster algorithm 450 adjusts the demodulation signals 
D\(k) to reduce or eliminate crosstalk. In one embodiment, the adjuster algorithm 450 
reduces crosstalk by configuring the demodulators, as discussed in more detail below. 

One skilled in the art will recognize that the lowpass filters provided in connection 
with the decimation blocks can provide other filter functions in addition to lowpass filtering. 
Thus, for example, the lowpass filters 420, 430, 440, and 450, and the decimators 420, 430, 
434, and 440 can provide other filter functions (in addition to lowpass filtering) such as, for 
example, bandpass filtering, bandstop filtering, etc. Moreover, the post-demodulation 
decimation rate R 2 need not be the same for each output channel. Thus, for example, in 
Figure 4, the decimator 440 can have a first decimation rate R 2 = n while the decimators 
430 and 434 have a second decimation rate R 2 = r 2 . 

The demodulators above are described in terms of digital signal processing on 
sampled data. Thus, the demodulator signals are written D\(k). The demodulators and the 
filtering associated with the demodulators can be done in using analog processing (using 
time-domain demodulator signals D\(t)) or on sampled data signals (using digital-domain 
demodulator signals D[(k)). For convenience, the following development describes the 
demodulator signals primarily in the time domain, with the understanding that the 
modulators can be implemented using digital signal processing or analog processing. 

The characteristics of the demodulation signals D\(t) and D 2 {t) affect how much 
crosstalk is seen in the output signals. In an diagonal system, that is, when the demodulator 

has been diagonalized, there is, ideally, no crosstalk. The first output signal S x (t) is an 

estimate (or approximation) to the signal S\(t). Similarly, the second output signal S 2 (t) is 
an estimate (or approximation) to the signal S 2 (t). When the composite signal S(t) is a linear 
combination of the signals S[(t), then the relationship between the signals S\(t) and the 

signals S g (t). When M x = cos cot , M 2 = sin cot , n(t) = Q, and there is no distortion (e.g., 
no pre-combination, combiner, or channel distortion) then: 

S(t) = S x (t) cos cot + S 2 (0 sin cot (3) 
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Then: 



If: 



then 



S } (t) = LP[D i (t)S(t)] (4) 

S 2 (t) = LP[D 2 (t)S(t)] (5) 

D^{t) = 2 cos cot (6) 

D 2 (t) = 2sincot (7) 



Z), (t)S(t) = 25, (0 cos 2 cot + 2S , 2 (0 sin cot cos ctrf 
= 5, (0 - £, (0 cos 2cot + S 2 (/) sin 2cot 



(8) 



After lowpass filtering to remove the terms with a frequency of 2 cot and higher 

S,(t) = S t (t) (9) 

Similarly for 

S 2 (t) = LP[D 2 (t)S(t)] (10) 

then 

D 2 (t)S(t) = 25, (0 sin cot cos art + 25" 2 (0 sin 2 cot 
= S 2 (r) - S 2 (0 cos 2cot + S i (t) sin loot 

After lowpass filtering to remove the terms with a frequency of 2 cot 

S 2 (t) = S 2 (t) (12) 

In the above analysis, it was assumed that there are no time delays or phase shifts in 
the signal S(t), and thus, configuration is relatively straightforward 

When an unknown delay (or phase error) is introduced, then the signals are no 
longer diagonal. Consider, for example, the situation when a delay A is introduced into the 
composite signal. Then: 

S(t) = cos co(t - A)5, (t - A) + sin co(t - A)S 2 {t - A) 

It then follows that: 
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5, (t) = LP[2 cos cot(cos co{t - A)S X (t- A) + sin co(t - A)S 2 (t - A))] 
= LP[2 cos cot(cos cot cos coA + sin cot sin (^A)^ (f - A)] 

+ LP [2 cos cot (sin cot cos <yA - cos cot sin <yA)5 2 (t - A)] 
= cos coAS x (t-A)- sin &>A5 2 (/ - A) 



The above equations can be expressed in matrix form as: 



"s,(0" 




cos<yA -sintyA 




A) 


J 2 (0. 




sintyA C0S6?A 




A) 



Then 



S^t-A) 
S 2 (t-A) 



coscyA sincyA 
-sinajA coscoA 



L 4(0. 



The above equation can be expressed as 





-A)" 




costyA sincyA 


■LP 


"A (05(0 


S 2 (< 


T A) 




-sincwA coscyA 


A (05(0 



= ZP 



cos <^A sin coA 
-sincoA coscoA 



A(0 



5(0 



A(0 

D a (0 



5(0 



where 
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Z), (0 = cos coAD { (0 + sin coAD 2 (t) 

= 2 cos coA cos + 2 sin &>A sin cot 

and similarly for Z> 2 (0 . Thus the modified demodulation functions D x (t) and D 2 (0 can be 

expressed as a linear combination of basis functions. If the time delay A can be predicted, 
then the demodulator functions can be calculated and programmed into the communication 
system. However, in many cases the time delay A is not known or changes over time. As 
described below, the demodulator functions can be determined by system calibration 
procedures. 

When an unknown phase shift (or phase error) is introduced, then there may be 
crosstalk in the system. Consider, for example, the situation when a phase error (j> x occurs in 
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the signal S\ (t) and a phase error fa occurs in the signal Si{i). The phase errors can be 
caused by intrinsic properties of the components, intrinsic properties of the system, 
component variations, time delays, etc. In the presence of the phase errors: 



A (t)S(t) = 2 A s \ (0 sin(tftf + fa ) sin cot + 2A 2 S 2 (t) cos(cot + fa ) sin cot 
= 24 Sj (/) sin tftf [sin cot cos ^ + cos cot sin 0, ] + 

2 A 2 S 2 (t) sin cot[cos cot cos fa +sinotf sin^ 2 ] (16) 
= A l S l (0 cos fa - A } S l (t) cos fa cos 2atf + A X S X {t) sin ^ sin 2cot 
+A 2 S 2 (t) cos fa sin 2#tf - (/) sin fa+A 2 S 2 (t) sin cos 2cot 



Thus, after lowpass filtering 

S } (t) = A x S x (t) cos fa - A 2 S 2 (t) sin ^ 2 



(17) 



The above equation shows crosstalk because S, (t) depends in part on components 
of S 2 (t) when A 2 * 0 and cj) 2 * nn where n = 0, ± 1, ± 2... . 
Similarly, 



S 2 (t) = ^(Osin^ + ^^(Ocos^ 
The above equations can be expressed in matrix form as: 



^cos^ -A 2 sin</> 2 
A x sin^, A 2 cos</> 2 J 



S 2 (0 



After inversion 



S 2 (t)j cos(^,-^ 2 ) 



cos <f> 2 sin$ 
sin fa 00s fa 



5,(0 
S 2 (t\ 



S 2 (t) 



1 



cos(^, -<p 2 ) 



cos ^ 2 sin^, 

sin fa cos fa 

A 2 A 2 



A(0' 
LAW. 



5(0 



(18) 



(19) 



(20) 



(21) 
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Then 



"5,(0" 




~A(')" 


S 2 (t) 




_A(0 



S(t) 

Where 5^ (0 are modified demodulation functions, given by: 

cos ^ 2 sin cot + sin^ cos cot 



(22) 



AC) 
AO j 



i 



cos(^ -&) 



- sin (j> 2 sin atf + cos ^, cos ctf 



(23) 



The modified demodulation functions have the form 



A(0 = 2>/M0 

7=1 



(24) 



Choosing the coefficients a } to eliminate crosstalk configures the modulator. 

In one embodiment, the coefficients a } can be fixed coefficients computed using 

known properties of the system. However, as system properties change over time, such fixed 
coefficients may lead to unacceptably high levels of crosstalk. Moreover, variations from 
device to device may cause the fixed coefficients to give unacceptably high levels of 
crosstalk. 

Higher performance (that is, lower crosstalk) can be obtained by computing the 
coefficients as part of a system calibration or initialization procedure. Such calibration can 
be performed during system startup (e.g., when the system is turned on, or when the system 
begins processing data, etc.) and/or at regular intervals. In one embodiment, the adjuster 
algorithm 450 computes the coefficients a } from the calibration data. In one embodiment, 

the coefficients a } are chosen by making four calibration-type measurements to measure 
four parameters , £ 12 > £ 2 i > m d £22 > where: 



I 

£n =S.|s,=as 2 =o =M COS0, 
£12 =5,L.o,s 2 =* = -Ba 2 sinfa 



(25) 
(26) 
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£>1 =S 2 \s l =A,S 2 =0= Aa i^ n< / > l 

£22 = S 2 s,=o,s 2 =fl = Ba 2 cos0 2 



(27) 
(28) 



where A and B are amplitudes. Then 

1_ A 

tan<*,=|^ 
tan 0 2 =-^ 12 



(29) 

(30) 
(31) 

(32) 



b22 

From the above equations, it is evident that crosstalk depends only on <f> x and (f> 2 . 
Moreover, ^ and ^ 2 can be chosen to eliminate crosstalk without knowing A or B. This is 
useful for systems such as pulse oximetry systems where absolute measurements of a 
channel are difficult or impractical, but where relative measurements (e.g., channel-to- 
channel measurements) are practical. 

The demodulation signals A(0 can be generated using the values of 0, and <p 2 from 
the above equations. Alternatively, the demodulation signals D x (i) can be generated from 
quadrature components as: 

A (0 = b\i si n A* + b l2 cos cot (33) 
D 2 (t) = b n sin cot + b 22 cos cot (34) 

where the coefficients b\j are computed from <j> x and <f> 2 . 

In one embodiment, the demodulation functions are adapted from baseline 
coefficients, which are then improved through a calibration or initialization procedure to 
produce actual coefficients. The baseline coefficients are typically obtained from known 
properties of the system. The actual coefficients are usually relatively close in value to the 
baseline coefficients. This provides one way to assess the operational status of the system 
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and to evaluate the calibration procedure. In one embodiment, if the actual coefficients are 
too different from the baseline parameters then it is assumed that the calibration procedure 
failed in some manner or that the equipment has failed in some manner, and appropriate 
measures can be taken (e.g., alert the operator, sound a warning, etc.) 

To find the actual coefficients, the demodulation functions are initially given by: 

D x (t) = a n sin cot + a l2 cos cot (35) 
D 2 (t) -a 2X s\ncot + a 22 cosoot (36) 

Where the coefficients a y are the baseline coefficients determined from known or 
assumed properties of the signal S(t). For example, in one embodiment a i} = 5 y . In one 
embodiment, where initial estimates are available for fa and (f> 2 , then the values of a i} can 
be computed as discussed above. 

The crosstalk reduction obtained using demodulation functions based on the 
coefficients a i} can often be improved by computing new coefficients a i} and 

corresponding new demodulation functions D.(t) where: 

D x (t) = a x ! sin cot + a n cos cot (37) 
D 2 (t) = a 2l sincot + a 22 coscot (38) 

The process of finding the coefficients a i} begins by measuring two data sets, x\(t) 
and X2(i), as follows: 

*i(t) = S(t)\ Si=AtS2=Q (39) 
* 2 (0 = S(0 W2 =* (40) 
The data sets x\(t) and x 2 (t) are used to enforce the following constraint: 

nT 

\x i (t)D i (t)d = 0 (41) 

o 
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where i = 1, 2, « = 1, 2, 3..., and 7 is a time period corresponding to one complete 
modulation cycle. From the above constraint and the definitions of the demodulation 
functions, it follows that: 



a u Jxj (t) sin cotdt + <2 12 (f) cos cotdt = 0 

0 0 
nT nT 

a n \x 2 (t)sincotdt + a 22 jx 2 (t) cos cotdt = 0 



It is convenient to define 



and to define 



where 



Then 



y u = Jx 1 (/)sin^^/ 

0 
0 

nT 

o 

nT 

Yu = jx 2 (t) cos cotdt 



r 9 



ni/2 



(42) 
(43) 



(44) 
(45) 
(46) 
(47) 



(48) 



(49) 



(50) 
(51) 
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In one embodiment, to reduce crosstalk, it is desired to find the coefficients a y 



closest (in the sense of minimizing some specified error, such as, for example, a least 



The term in parentheses can be described as the baseline crosstalk. 

One of ordinary skill in the art will recognize that optimization methods other than 
least squares can be used. The solution methods for configuration are, for simplicity, 
described above in terms of a two-channel system. Using the above teachings, the extension 
to multi-channel systems is straightforward. 

Although described above in connection with a particular embodiment of the present 
invention, it should be understood the description of the embodiment is illustrative of the 
invention and are not intended to be limiting. Although described above in connection with 
a pulse oximetry system wherein a parameter to be measured is the attenuation of red and 
infrared light passing through a portion of a subject's body, it should be understood that the 
method and apparatus described herein can also be used for other measurements where two 
or more signals are passed through a system to be analyzed. In particular, the present 
invention can be used to demodulate two combined parametric signals responsive to the 
system to be analyzed where the two parametric signals have a predetermined timing 
relationship between them, as described herein. The invention can be used in connection 
with various physiological parameter measurement systems, such as, for example, systems 
that measure blood constituents, blood oxygen carboxyhemoglobin, methemoglobin, 
glucose, etc. Various modifications and applications may occur to those skilled in the art 
without departing from the true spirit and scope of the invention as defined in the appended 
claims. 



squared error) to the coefficients a i} such that the above constraints are satisfied. One 
solution, obtained by minimizing the least squared error is: 




(52) 
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